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Abstract 

The conformations of interacting linear polymers on a dynamical planar random lattice 
are studied using a random two-matrix model. An exact expression for the partition function 
of self-avoiding chains subject to attractive contact interactions of relative strength 1/c, < 
c < 1, is found as a function of chain length L and c. The number of configurations L a e bL 
as L — > oo is determined, showing that a chain undergoes a third-order collapse transition at 
c = y/2 — 1; the universal scaling laws are found. 
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1 Introduction 



A long polystyrene chain molecule in a solvent undergoes a coil-to-globule transition M at a 
finite temperature T = due to attractive polymer-polymer interactions. Many approximate 
theoretical techniques have been developed to describe this and other processes of macro- 
molecular folding Pj. In this letter exact results for this problem of polymers, modeled by 
self-avoiding chains with contact interactions, will be given on the ensemble of two-dimensional 
simplicial lattices with random fluctuations of the local intrinsic curvature which appear in the 
study of two-dimensional quantum gravity [||, |j] . This ensemble represents an annealed average 
over a certain class of fractals. The advantage is that statistical mechanics problems are often 
exactly solvable in this case using random matrix models which yield expressions for the 
thermodynamic and some of the geometrical quantities of interest. Moreover it is intuitively 
plausible, and found to be the case in all of the many solved examples, that the phase structure 
of such statistical systems is the same as on a fixed lattice but with weaker transitions when 
one averages over fluctuating lattices. Although a simple structureless polymer will be studied 
here, the same techniques could be extended to more complicated polymers types. 

A two-dimensional simplicial lattice is obtained by gluing, pairwise along edges, A squares 
each of edge length 6, so as to form a surface with spherical topology say. A polymer chain 
of L links walks isotropically on the edges of squares, cannot cross itself (excluded volume), 
but can occupy a given edge any number of times (figure 1); the latter represent contacts. 
The partition function is a sum over all possible simplicial lattices and configurations of a 
polymer on such lattices with a weight for contacts. The G-temperature is found exactly in 
this model and represents a third-order collapse transition. One is therefore led to predict 
that on a fixed two-dimensional lattice (e.g. honeycomb), with the rules of contact given in 
Fig.l, the transition will be first or second order.[] The universal exponent a, defined by the 
number of configurations as L — > oo, L a e bL , changes from a = —5/2 at T > to a = —3/2 at 

1 Some previous exact results p| on fixed lattices have used a different definition of contact, forbidding entirely 
configurations 'a' and 'b' of Figure 1; a contact occurs when two non-adjacent parts of the polymer come within one 
lattice spacing. This could make a difference in the collapsed regime since the polymer can be much more compact 
if multiple occupation of an edge is allowed (type 'a' in Figure 1). Such multiple occupation necessarily occurs, for 
example, in the high-temperature expansion of the lattice 0(n) model with action S.S 
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T = 0. For T < it was not possible to determine a. The universal behaviour around T = 
is described in terms of a reduced temperature T — ~ L -1 / 6 . 



2 Matrix Model. 

The self-avoiding walk on fluctuating random surfaces was described using a one-matrix model 
by Duplantier and Kostov ||. In order to include an attractive interaction whenever contacts 
occur consider a two-matrix model with hermitian NxN matrices (ft and tft, 



VcftVipexp yNTr 



1 j.2 i 1 ,2 / / 5(1 C 2 ) 2 4 2\L-n / 2L 

2<A + 2^ ~ c< ^ 4 ^ + I 1 - c ) r ^ 



(1) 



As usual, the Feynman diagram representation of the expansion in g generates graphs dual 
to simplicial lattices and as N — > 00 those of genus zero dominate. T is a chemical potential 
conjugate to the number of polymers n, each of length L, i.e the surfaces are built from (ft 
vertices dual to squares while the polymers from tp matrices. The last term in in the action 
TrV((ft, ip) represents a hole of length 2L in the surface with one point on the boundary marked; 
one can visualize the polymer by sewing up the hole starting at the marked point. Each possible 
hole and each possible location of the marked point give a different polymer configuration. The 
polymer is connected into the surfaces by (fttft so that c represents a polymer-lattice Boltzmann 
weight, while the polymer-polymer (contact) weight tft 2 has been normalised to one. Thus 
< c < 1 sets the temperature scale, with c = 1 being infinite temperature. 

Z N ^oc= £ c w T n g A H(A,w,n,L) (2) 

w,n,A 

where the sum is over genus zero surfaces and H is the number of configurations with n 
polymers, each of length L with w polymer-contacts, on all surfaces of area A. Since H grows 
as ~ const" 4 there exists a critical g = g c at which large area surfaces become significant and 
a universal continuum limit can be established Q. 

The matrix model ([!]) can be solved in the large N limit by the method of orthogonal 
polynomials 0. The polynomials Pi((ft) = (ft 1 + • • • , Qj{ift) = 1ft 3 + ■ ■ ■ are such that 

»+oo 

dcftdifte-V^Pii&Qj^) = hidij (3) 



and the free energy is — logZ = const — ^^^(N ~~ n)log fi, where fi = hi/hi-i and the 
constant is ^-independent. P and Q satisfy recursion relations 

cpP,{4>) = P i+ i + r i P i - 1 + s i P i „ 3 + --- (4) 

^QjW = Qj ■ I • <l.jQ.i I • l,Qj 3 (5) 

That for P terminates after L + 1 terms. From identities such as 

##e _y( ^ ) ^-iQi = (6) 



one can find a closed set of equations for r, s, q, t, and / which can be solved for /. In the large 
N limit one can define a continuous variable x = i/N and functions fi = Nf(x),ri = Nr(x), 
etc. and these equations are explicitly: 

cr = f-Ty/f6[i,i-l] (7) 

cq = f-3g(l-c 2 ) 2 fr (8) 

ct = -g(l -c 2 ) 2 / 3 (9) 

x = -cf + q-T^e[i-l,i] , (10) 

where 

0[i,j] = 2L(1 - c 2 ) L I d0dV>e-^)-p=^ 2i - 1 -%= . (11) 

J \hi ^Jhj 

Substituting eqs. (0) and (|) into eq. @ gives an equation (the 'string equation' in quantum 
gravity) for /. will be evaluated shortly for large N. 

The critical point for continuum surfaces is identified by setting T = (no polymers) and 
finding the singularity in the solution for /(l) of the quadratic equation ( |i"0| ) as g is varied, 
yielding g c = and f c = 2c/ (1 — c 2 ). Defining <7 = g c (l — fid 2 ), where 5 is the lattice 
spacing, \i is the renormalised 2D 'cosmological constant' conjugate to renormalised surface 
area A = A5 2 . In the limit 5 — > 0, \i couples to universal continuum surfaces (A finite) and 
surfaces finite in lattice units (.4 = 0). The partition function will have a part regular in /x, 
usually representing the latter surfaces, and a non-analytic part representing the former. It is 
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appropriate to define scaling variables x = 1 — z5 2 and f(x) = f c (l — u(z)8) in which case, in 
the limit 6 — ► 0, 

u 2 = fi + z (12) 



from eq,(|lC)|) at T = and 

-\ogZ = N 2 5 5 / dz zu + regular (13) 
J o 

3 Polymer Collapse. 



First the case of a single polymer will be addressed, when eq.(12) is appropriate and one 
studies the expectation value of Tr(t/; 2L ) in this theory, i.e. Z expanded to first order in T, 

Z = Z + TZ 1 + --; 



Zi = [v<f>Vif>(l - c 2 ) L Tr(^ L )e- 7VTr ( v '( r = )) 

= (i _ jfy / ^-vtN^ift . ( M ) 



Thus, if Z\ = Y^i=i Z'xi from eqs.(||)(||)(||) and (Q) one finds 

= ^.^fe^a^ (15) 

(^-2p)!(p + L)!p! 

^ = i-^(l-^ + 0(5 2 )) 

C D 

([L/2] indicates greatest integer < L/2). In order that the polymer be an extended object on 
the surfaces of finite A one must scale its length as L = 1/5 in which case 

Z[ = G(L, c ) e -Mc)ui+0(S) + ( CO nst~ ^) (16) 

where G oc L a a b ^ L . Despite many attempts the author has not managed to derive the asymp- 
totic form ( |l6|) for all c. It is generally expected however that the number of configurations 
has the form of G as L — > oo, which numerical evaluation of (|l^) confirms. Figure 2 plots A(c) 
obtained from 

Z[(L,c,0) 



ulX = log 



Z[(L, c, ul) 



(17) 



for large but finite L — corrections to the right hand side above are 0(5). 

G = linif^oo Z[(L, c, 0) is a kind of wavefunction renormalisation which represents the 
effects of surfaces of zero renormalised area, A ~ 0{5 2 ); it does not depend upon u or fi. On 
the other hand u\ represents an induced free energy per unit length of the polymer due to 
fluctuations of the finite A surfaces; it is non-analytic in [i from eq.(|l^). The behaviour of 
A(c) thus governs the interesting geometrical properties of the polymer. In particular, fig. 2 
seems to show a transition to A(c) = at c ~ 0.4, which will be verified shortly. Physically, 
as c is reduced the polymer favours polymer-polymer rather than polymer-solvent bonds and 
so the induced free energy A falls. Below a certain c c ~ 0.4 the polymer collapses entirely and 
is only connected into the surface by a small number of polymer-solvent bonds. This is the 
coil-to-globule transition on fluctuating random surfaces. The polymer-polymer contacts also 
induce free energy per unit length but at a lower order in S, and it will be necessary to choose 
a different scaling law for L in order to see this contribution near c = c c . 

In order to find c c exactly, determine the universal exponent a, and the details of the phase 
transition, it is useful to rewrite Z[ as 



Z[ = (1 - c 



2\L 



n dp 

-7T 2-7T 



V?e- ip + -4^e ip + 



t 



'Sip 



2L 



(18) 



v7 P' 2 

As the continuum limit is approached it is expedient to rescale p — ► p6 d for some exponent d. 
Expanding in 5 one has 

ip = — [a x + a 2 ip5 d + a 3 u5 + a 4 p 2 5 2d + a 5 iup5 d+1 + a 6 u 2 5 2 + a 7 fi5 2 H ] (19) 

Consider first the high temperature c — > 1 limit. It is known from the c = 1 solution that 
d = 1/2, so this ansatz will be used in the high temperature (coil) phase. Writing c = 1 — e 
for small e, the coefficients en can be found as power series' in e; in particular a 2 = 0(e) and 
one finds from fllg|) 



Z[ = V 4^ [8 + ° (e)]LeXP {~ UI i 1 " I + ° (e2) )) (20) 
in the continuum limit 5 — > 0. This is of the form (16) with a A that correctly matches onto 
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the answer of Fig. 2. The partition function itself is then 

ZiaA ^=^+i.) e -A<V* (21) 

using eq. (|l2|) . Large area A — > oo corresponds to [i — ► 0, in which case the universal power 
exponent is a = —5/2. If c is reduced by a finite amount from one, 02 is finite and the 
resultant oscillatory behaviour means that is no longer possible to neglect the higher orders in 
p as 5 — > 0. Nevertheless one expects the form (^Tj) to persist down to c = c c . 

Fortunately there is one other value of c where 02 vanishes, given by c = y/2 — 1; 03 also 
vanishes at this point indicating that this is the collapse temperature. A non-trivial scaling 
solution is obtained with the scaling laws c = y/2 — 1 + crS 1 / 3 , L = l/d 2 , and d = 1 (note that 
polymer length scales as area now). In this case one derives from fll8| ) 

plfJ. poo 

Z x oc const. L N8 z - r / dy e -^-^ 3 Vvl (22) 

in the continuum limit, where a ~ 41; a = —3/2 at a = 0. As the reduced temperature 
a — ► +00 one recovers the form (21) with A — > acr 3 matching in the high temperature phase 



(fig. 2). At g = 0, Zi, while now analytic in fj,, in this case seems to represent finite A surfaces. 
A = surfaces — A ~ 0(5 2 ) more precisely — are characterised by the fact that a finite 
number of derivatives with respect to fi can remove their contribution, since each derivative 
brings down a factor of A. But Z\{a = cannot be removed by differentiating. It can be 
interpreted as a surface dense with polymer whose area is proportional to the polymer length; 
for small I differentiating does remove its contribution. It is tempting to extrapolate into the 
collapsed phase a — > — 00 also, but the above scaling laws do not seem to be consistent since 
Z\ becomes a saddle point integral with the result Z\ ~ (u 3 /l)exp(fil). This cannot be correct 
because Z\ should be monotonically decreasing with fi. 

If polymers are present with fugacity T one must go back to eq.(|i~2"D and modify it accord- 
ingly. From eq.(||) and eq.([To|) one has 



x = -cf + L - ^ + 69f ^ T e[i, i - 1] - 2 J]T6[i - 1, i] (23) 

In order that the new terms contribute to the scaling part of the equation, but do not disturb 
the critical condition for continuum surfaces, g c = i to leading order in 5, one must tune 



r — > appropriately. For c > y/2 — 1 one finds at order S 2 in eq. ( |23| ) 

= ^ + z-n 2 - -f\V~le- lXu . (24) 

The wavefunction renormalisation G of Zi has been absorbed in T and a renormalised fugacity 
7 defined, of dimension <5 5 / 2 . When 7 ^ the singularity at g = g c (/i = 0) is now shifted 
by an 0(5 2 ) term, \i = \i c say, which therefore couples to renormalised area A. \i c represents 
the free energy of the polymer system per unit A on finite surfaces. The positive root of 



eq.(24) is the correct one since u — > -\-y/fl + z as 7 — > 0. A singularity in u appears when 
2u = r y\ 2 l 3 / 2 e~ lXu . Then u(p) is singular at [i = /i c . In particular at c = y/2 — 1 + e, e > 0, 
one finds u(/j c ) = ±a 2 ~/l 3 / 2 e 6 + 0(e 7 ) and y, c = ajVk 3 + 0(e 4 ). For e < 0, A = = /x c , 
so that 3 3 /i c /9c 3 is discontinuous implying that the collapse transition is third order The 
numerical factor a is not expected to be universal but the 7 and I dependence should be. 
Finally, it would be interesting to study the possibility of inter-polymer aggregation in the low 
temperature phase. 
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FIGURE CAPTIONS 

Figure 1. — A portion of a two-dimensional simplicial lattice made of square simplices (broken 
lines) with the polymer chain shown as a continuous line. This figure shows one contact 
interaction on edge 'a' ('b' is not counted as a contact). 
Figure 2. - The induced free energy per unit length u\ estimated from eqs. 

©(0) at L = 150. 

Broken lines are the leading perturbative results around c = 1 and c = y/2 — 1 for L = oo. 
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